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STRESSES IN A TWO-BAY NONCIRCULAR CYLINDER UNDER TRANSVERSE LOADS 


By GEORGE E. GRIFFITH 


SUMMARY 


wi method, taking into account the effects of flexibility and 
based on a general eighth-order differential equation, ts pre- 
sented for finding the stresses in a two-bay, noncircular cylinder 
the cross sectiun of which can be composed of circular arcs. 
Numerical eramples are giren for two cases of ring flexibility 
for a cylinder of doubly symmetrical (essentially elliptic) cross 
section, subjected to concentrated radial, moment, and tan- 
gential loads. The results parallel those already obtained for 
shells with circular rings. 


INTRODUCTION 


In airplane fuselages with flexible rings subjected to con- 
centrated loads, the stresses in the neighborhood of the load 
differ markedly from those given by the simple engineering 
formulas, and more refined methods, which take into account 
the interaction of rings and shell, are needed to predict the 
stresses accurately. The first paper on this subject, pub- 
lished in May 1944, was that of Wignot, Combs, and Ensrud 
(reference 1}, who treated the circular cylinder subjected to 
concentrated loads but neglected the effect of the extensional 
deformations of the shell. Hoff (reference 2) gave a more 
complete analysis, including the effects of many rings, for 
the case of symmetric transverse loads. The results were 
corroborated experimentally by Kuhn, Duberg, and Griffith 
(reference 3), who also extended the theory to include con- 
centrated moment and tangential loads. Later, Duberg and 
Kempner (references 4 and 5) reduced the labor of computa- 
tion by giving the results in the form of charts and showed 
that for practical purposes it was usually sufficient to con- 
sider only a region within 2 bay lengths of the load. Further 
investigations considered additional effects, heretofore neg- 
lected, such as the shearing and axial deformation of the rings 
(references 6 to 8), bending rigidity of stringers (references 
8 and 9), shear carried by the stringers (reference 7), and 
eccentricity of ring and sheet (references 6, 7, 8, and 10). 

All investigations referred to dealt exclusively with rein- 
forced monocoque cylinders of circular cross section. The 
present report gives an analysis for a two-bay noncirculer 
cylinder, enclosed between ring bulkheads rigid in their 
planes, with the middle, flexible ring subjected to concen- 
trated and distributed loads. The fundamental assumptions 
used in the analysis are the same as those previously used for 
circular cylinders. 


Many noncircular fuselages can be closely approximated 
by using circular sections of different radii and joming the 
sections at points of tangency. The rings discussed herein 
are of this form. Associated with each ring section is a two- 
bay panel (fig. 1), any number of similar panels composing 
the complete two-bay cylinder. The solution for the stresses 
in such a structure is based on the development of a general 
eighth-order differential equation, written in terms of the 
moment at the skin center line. A separate differential 
equation of the same general form applies to each curved 
panel in the structure. Application of the correct boundery 
conditions results in sets of simultaneous equations which 
yield the unknown constants in the moment expression. 

Although a two-bay cylinder does not conform to the usual 
fuselage structure, it is believed, on the basis of comparisons 
with some of the work previously mentioned, that the results 
obtained are indicative of those found in more complicated 
structures. In accordance with the findings of prior mvesti- 
gations, shear and axial deformations of the loaded ring are — 
neglected, but eccentricity of ring and sheet is included. 

The numerical examples deal with a doubly symmetrical, 
two-bay cylinder of nearly elliptic cross section (fig. 2), sub- 
jected to concentrated loads. 


SYMBOLS 


Ci coefficient of shear flow in sheet 

Ca coefñcient of ring axial force 

Cy coefficient of ring bending moment 

Cr coefficient of ring shear force 

E Young’s modulus in sheet and ring, pounds per 
square inch 


„Rigid 





FIGURE 1.—Two-bay panel showing positive directions of displacements and sheet stresses. 


1 Supersedes NACA TN 2512, “Stresses in s Two-Bay Noncircular Cylinder Under Transverse Loads” by George E. Griffith, 1951. 
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FIGURE 2.~—Two-bay s deart wa composed of four circular ri poet ll with doubly symmetrical 


G 
iT 


anelreular rings, used for numeri 


shear modulus in sheet, pounds per square inch 

axial force acting on sheet, pounds 

moment of inertia of ring section, inches ‘ 

length of bay, inches 

bending moment at sheet center line, inch-pounds 

applied concentrated moment, inch-pounds 

bending moment about neutral axis of ring, inch- 
pounds 

applied concentrated radial load, pounds 

radius to sheet center line of circular panel, inches 

radii of circular panels composing cylinder, inches 

applied concentrated tungential load, pounds 

shear force, pounds 

complex roots of auxiliary algebraic equation (used 
when real roots also occur) 

complex roots of auxiliary algebraic equation (used 
when no real roots occur) 

eccentricity of ring and sheet (distance between 
sheet center line and neutral axis of ring), Inches 

distributed radial load acting on sheet, pounds per 
inch 

distributed axial load acting on sheet, pounds per 
inch 


real roots of auxiliary algebraic equation 

distributed moment acting on sheet, inch-pounds 
per inch 

normal force in stringers, pounds per inch 

shear flow in sheet, pounds per inch 

distance in circumferential direction, inches 

thickness of sheet, inches 

thickness of all material carrying bending stresses 
in panel if uniformly spread around circum- 
ference, inches 

displacement of sheet in axial direction, inches 


v displacement of sheet in circumferential direction, 
inches 
to . displacement of sheet in E direction 
_ atring, inches 
wW displacement of sheet or ring in radial direction, 
“ inches 
T distance in axial direction measured from baia 
ring, inches 
y shear strain in sheet 
P normal strain in sheet in circumferential direction 
Ex normal strain in sheet in z-direction 
6 angular distance, degrees or radians 
g : stringer normal stress, pounds per square inch 
T: “normal stress in sheet in circumferential direetion, 
` pounds per square inch 
T ~ -shear stress in sheet, pounds per square inch 


GENERAL DIFFERENTIAL EQUATION 


The basic element of the present analysis is a two-bay 
panel, as shown in figure 1, composed of sheet, longitudinal 
stiffeners or stringers, and transverse stiffening ring sections 
of constant radius. Distributed or concentrated loads are 
applied in the plane of the middle or flexible ring. By 
joining several panels at points of tangency, many two-bay 
cylinders of various cross-sectional shapes can be achieved. 

In actual practice the sheet covering is outside the rings, 
and this eccentricity of sheet and ring is henceforth taken 
into account. For convenience the moment considered is 
the moment which exists at the sheet center line rather than 
at the neutral axis of the ring. From this moment the 
bending moment in the ring is easily found. 

The underlying assumptions used in the analysis are as 
follows: 

(1) The shear stress, carried by the shect alone, may vary 
in the circumferential direction but remains constant in the 
axial direction. 

(2) The material in the cross section of the panel (sheet 
and stringers) capable of carrying normal stresses duc to 
bending of the panel is assumed spread around the cireum- 
ference in a fictitious sheet of thickness t’. 

(3) The loaded ring has no torsional stiffness or bending 
stiffness out of its plane. 

(4) The end ring supports are restrained from deforming 
in their planes but are free to warp out of their planes. 

Under these assumptions, for any panel with constant 
geometrical properties, a general differential equation is 
developed for the moment at the sheet center line. All 
forces, stresses, and displacements in the panel may be 
obtained from that moment. The general differential equa- 
tion (see the appendix for developinent) is 


4A, 
ete 3B) opt (168+ ee oe 


a| BEE (2—F4 a) og + 6A(1— ati) M= F(6) 
(1) 
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where 


FO=R| pt0- 3B) Se 13S + 
f 2a) df 
e| g+- 3B)54—-3(B—75 cf 


(i ate) |+ 
Se z d?m 


Rete 3D) gr r+- Baa dë 


a| B-30 +5) a 


When the ratio e/R, a measure of the eccentricity of the 


ring and sheet, is very small ER 3) it ean be neglected; 


equation (1) then contains only the two nondimensional 
parameters <1 and B, which are related to the geometry and 
PR? Et R? 
A= 77s and B= Gre 
The ratio A/B is an index of the over-all flexibility of the 
structure. When this ratio is very small the ring section is 
rather rigid. A large value of A/R (500, for example) indi- 
cates that the ring section is somewhat flexible, so thet radical 
departures of the stress distributions from the elementary 
values result. In actual practice the parameter B varies 
over a small range (between about 10 and 80), whereas A 
is usually much larger and may be as great as 2X10", or 
even more; consequently, flexibility may be thought of in 
terms of parameter .i alone. Values of A less than about 
200 usually indicate relatively inflexible rings, and increasing 
values indicate increasing ring flexibility. 
The right-hand side of equation (1), F(6}, contains terms 
resulting from the application of distributed loads; if only 
concentrated loads are present F(#) becomes zero. 


physical properties of the panel: 


S13 
SOLUTION OF GENERAL DIFFERENTIAL EQUATION 


The solution of equation (1) depends upon the nature of 
the roots of the auxiliary algebraic equation associated with 
the differential equation. Of the required eight roots of the 
algebraic equation, four may be real and four complex, or 
all eight may be complex. (For the special case of no 
eccentricity of ring and sheet, sets of roots for several com- 
binations of parameters A and B are given in table I.) 

If there are four real roots +k, +h, and four complex 
roots +(a+62), the solution for the bending moment is 


AL = Citt Oue 1? + Coe"? + Ciye 2? + 
Cre Ste Ope -0E Oe et 
Cyme ~ “~*? + Particular solution (2) 


which can be written in real form as 


=C sinh &,6+-C, cosh £,0+C; sinh 44+ C1 cosh k+ 
C, sinh af cos b8 +C, sinh a8 sin b8 +C; cosh ab sin b04 
C; cosh a8 cos 66+ 2-7 (6) (3) 


If there are four pairs of complex roots +(¢,+¢,t) and 
+ (cg+-dt}, the bending moment is 


M= Cye tiD L Cpe TAD One rte et CF ge Oe 
Cre este tL Oye 2 GOI Ope tD Cg eae. 
Particular solution (4) 


which, written in real form, becomes 


Af=C, sinh ¢,@ cos dif + C, sinh c0 sin d,¢+ C cosh c, sin d,6-+ 
C, cosh ¢,6 cos di0 + O; sinh e208 cos d,6+ C, sinh c8 sin d.6+ 
C, cosh c8 sin d+ C; cosh ¢,4 cos da0 +3A (0) (5) 


Since the particular solutions gJ(@) and 7/(@) depend upon 
the form of F(6}, no formal solutions are given here. Further 


TABLE I.—ROOTS OF AUXILIARY ALGEBRAIC EQUATION WHEN e=0 FOR SEVERAL COMBINATIONS OF A AND B 
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+ i. 250205 +0. 956180 | +0. 743842 +0. 571042 +k; 
SEET” +7. 745793 | 413. 416404 | +21, 213203 432. 403439 + ks 
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+2. 523294 
+ 13. 416028 
+ (L 236394+ 2. 3889510) 


+3. 157270 
10 +7. 728363 
+ (1. 447531 +2. 8183909 | 





+5. 016203 | +3. 835266 
10 +7. 539604 l +13. 412609 
4 (2. 058771 + 4 030199) | + (L 837706 +3. 4177220 
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+2 117014 + 1. 795392 +k, 
+21. 213188 +32. £03439 +k; 
+ (L 052130 +2. 0850722) + (0. 896469 + 1. 8457882) = (a + be) 
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+3. 230542 ) ! 
4-21. 213049 | 
+ (1. 59501742. 9595622) | 
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remarks concerning the solution of equation (1) are confined | 


to the solution in real form, either equation (3) or equation 
(5). 

Equation (3) or equation (5) expresses the moment in any 
particular panel where A, B, and R are constant; such a 
solution exists for each panel constituting the structure to 
be analyzed. For the case of a circular cylinder, solution 
(3) or (5) is the same as the energy solutions of references 
2 and 3 for a similar two-bey structure. 

RELATION BETWEEN THE MOMENT M AND OTHER FORCES 
AND DISPLACEMENTS 

All forces and displacements in the panel may be expressed 
in terms of the moment given by equation (3) or (5) and 
derivatives of this moment. As previously noted, this 
moment is the moment at the sheet center line, and with the 
exception of the ring bending moment all other forces and 
displacements are those at the center line of the sheet. When 
no eccentricity is involved the sheet center line coincides with 
the neutral axis of the ring. Of major interest in a structure 
such as the one described herein are the forces, stresses 
(obtainable from the forces), and displacements listed in 
the following paragraphs, together with their mathematical 
expressions which are readily obtained in general form in the 
appendix. These expressions become considerably simpli- 
fied in the absence of distributed loads. (For positive sign 
convention of the displacements and forces see figs. 1 and 
3, respectively.) 


f 


Am af 
Fh dð Á 
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Figure 3,—Free-body diagram of ring section showing positive directions of forces. 
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The eight quantities listed in this paragraph are associated 
with the boundary conditions (discussed in the next section). 
These quantities—the forces and displacements at the 
flexible ring, the shear flow at the panel edge, and the axial 
displacement between the flexible and rigid rings at the panel 
edge—are: 

(1) Bending moment, Af 

(2) Shear force 
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(3) Axial force 


1 ËM 
mp Gg FRI GG 


(4) Sheet shear flow 
1 (d@At, BAL d'm 
; -e| itg t py (m+ call 


1-373) do de 





(5) Axial displacement 


_ LL (ËM dM d?f dm d'm 
UTR de do ata (a tae) | 
(6) Tangential displacement 
Mf EM 


-r| r+ +7(m + oR) + 
Phys. za „d'm dn) 


g ——srai Get a “de 
D? sees M oa _ 
6 EV TOG TR 


(7) Radial displacement 





Z [EM UM pdh, df 1 (dm, dm 
w= SF e ae TAr -el t ata +o) |}+ 
L3 dM dM a dh d'm . dm 
meala Hae el S45 iti et ae) h+ 
e d AL dm 
pr (tp Ge BS 6 Ge 


(8) Rotational displacements 





dw vo  L (dM, ,@M, BM dth 
det aa t ae + ge | he G+ 





THR (m+2 Tat ++ sREP Ir + 


M EM pf Ph d'h, df, dy 
2 tage E Fear tat t+ 


gdm dim dif, eM df „Em 
2ga ae) her Ge ia TR ae EI ae 
The shear and normal stresses in the two supporting panels 
follow. The shear stress r in the sheet is given by 


maa tae” | ate (tae) $ 


The stringer normal stress ¢ is 


ËM , d'M df dm d'm 
= Sha ae qe t dg setae a aet de a) | i 
The ring bending moment—that is, the bending moment 


about the neutral axis of the ring—is 


Me=M+5O" enfe @ (8) 


For no eccentricity the moment at the sheet center line be- 
comes the ring bending moment. 
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BOUNDARY CONDITIONS 


Inasmuch as a general differential equation of the form of 
equation (1) applies for each of the panels composing «4 
given structure, a like number of bending-moment solutions 
of the form of equation (3) or (5) results. For the deter- 
mination of the unknown coefficients appearing in these 
solutions s sufficient number of boundary conditions must 
be found, one for each coefficient. Regardless of the number 
of panels, the boundary conditions involve only the eight 
expressions listed in the previous section, consisting of the 
forces and displacements which occur at the sheet center 
line. 

Although, for a particular problem, the application of 
these expressions depends upon the structure to be analyzed 
and the loading, all the boundary conditions may be 
summarized in one general statement: No change in 
displacement or forces can occur across a boundary unless & 
concentrated force is applied at the boundary, in which case 
the change in forces must equal the applied force. Any 
concentrated load in the plane of the ring can be resolved into 
a racial or shearing force, a tangential or normal force, and 
a moment. Then when a concentrated load is applied, the 
boundary conditions require that the difference in shear forces 
of the adjacent panels be equivalent to the applied radial 
load, the difference in axial forces be equivalent to the 
applied tangential load, and the difference in moments equal 
the applied moment. In the absence of any concentrated 
loads, all the forces and displacements must be continuous; 
that is, all eight expressions in one panel must equal the 
corresponding expressions in the adjacent panel. If a con- 
eentrated load is applied within a panel, it is necessary in 
the analysis to consider the point of application of the lead 
as a boundary and, bence, to consider the panel as two 
panels, one on either side of the load. 

Although terms associated with distributed loads appear 
in the expressions for the boundary conditions, distributed 
loads affect the Loundary conditions only indirectly inas- 
much as tbey affect the displacements. 

Further discussion of boundary conditions, as they apply 
to one of the numerical examples, is given in the following 
section. 


NUMERICAL EXAMPLES 


The numerical examples deal with two geometrically 
similar cylinders constructed of four panels forming a 
doubly symmetrical, essentially elliptic cross section (fig. 2), 
with the flexible ring subjected to concentrated radial, 
moment, and tangential loads at an intersection with the 
major axis. Cylinder 1 has & very stiff loaded ring and 
cylinder 2, a relatively flexible loaded ring. For each cylinder 
the sheet thickness is constant, there is no eccentricity of 
ring and sheet, and the radius R, of the top and bottom 
panels is one-third the radius R+ of the middle panels. The 
moment of inertia of each ring is constant, but because of 
the change in radius, the rings change in relative stiffness 
from one section to another as indicated by the change in 
1/B given in table II. As seen in figure 4, the top panel, 


TABLE II.—_NUMERICAL VALUES USED IN COMPUTATIONS 























URS „EYR 
| 4= 75: B= HIE 
PR - | A Roots of auxiliary 
der | section A B 5 equation 
| | | | | +0.238301 
--50 | 20 
land3; zo lF + 2.581786 
| i + (0.109747 + 1.019532 
i 
+1.061451 
2 50 | 20 +7.745880 
+(0.521510 + 1.354493 
| 150000 | 20 +-(2.509792 + 0.621665) 
Land 3 i=} 925 
pees 729 | 9 ° | 4.(0.915513-4 2.101999) 
2 z 
. 4+ 5.534415 
2 150,000 | 20 | 7500.00 + 7.395633 
| (2.181206 + 4.291 882i) 
@=0° 





20° 


S 
Section 3 


| 
1go* 
FIGURE 4.—Cross-sectional view of half of loaded ring used for numerical examples, 


section 1, joins the middle panel, section 2, at @=60°, and 
the middle panel joins the bottom panel, section 3, at 
§=120°. (These dimensions were also used in constructing 
fig. 2.) 

Comparisons of the calculated distributions of bending 
moment, shear force, and axial force in the ring and of 
shear flow in the sheet with the distributions given by the 
elementary theory are shown in figures 5 to 10. The neces- 
sary numerical values used in the calculations are given m 
tables II and III. For ease in reading figures 5 to 10, the 
abscissa, although it actually represents distance slong the 
perimeter, is given in degrees measured from the vertical 
axis of symmetry (as shown in fig. 4). Thus, since the ring 
perimeter of section 2 (fig. 4) is three times that of section 1 
or section 3, whereas the angular distance for all three 
sections is the same (60°), the distance in figures 5 to 10 
along the abscissa from 120° to 60° (corresponding to 
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TABLE III.—VALUES OF COEFFICIENTS! FOR NUMERICAL EXAMPLES | 





Cylinder 1 


Coefficient 





‘Load 
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_ Cylinder 2 














































Moment Tangential Moment Tangential 

See Sees Seer . 0. 147507 0. 249767 i rele aaa —O. 860733 —03 0. 119301 —03 
Daaa- |. 27504104 | 15740603 | we |. 29273103 | — 371518—04 
seeccnenatcenns | 418616 =| lk a | asemaa | 14801801 | —. 22123802 
MNE ear O0ONS B © 286338—03 Daanan | — 521539—02 | 481959 —03 
ecommerce es , 354900 —02 EE Gece ee ee eer for ree 
. 171565 —01 . 188166 . 202467 —. 657213 —08 . 140829 —06 — 161423 —07 
Cio . 303490 —01 —. 121785 —, 111398 l —. 534558 —01 . 270768 — 697692 —01 
Cy —. 264559 —10 . 116591—10 —.116591—10 |, ._ .397640—09 —. 161540 —08 . 710774—09 

Cie . 889303 —. 508263 —01. . 851830 . —: 479259 — 380612+ 01 . 638771 
Cu . 329176 —01 — 296012 . 472934 —. 122989 —O01 —: 361474—01 . 485865 —02 
Cu —. 490646 - 373871 —01 . 290987 . 829422 —03 — 491296 —O1 . 702908 — 02 
Cis . 192013 —. 184225 — 160403 — 117551—02 . 471184—O1 —, 668234—02 
. 720834—O1 - 274929 —. 400188—O1l 7 - 131739 —01 . 392728 —01 —. 527300 —02 
—. 119460 —. 164389—O1 , 174960 | =, 119054—05 - 778086 —0 £ — 192398 —04 

Cis . 315057 . 363960 . 224540 . go 202194401 —. 711069+ 01 . 872951 
Cio —. 898224 —08 . 174922 —09 . 301271 —07 . 342801 —04 . 801877 —04 . 168139 —04 
Crp —. 711275 —02 —. 259599 —. 364881 . 216297+01 . 882602 +01 —. 1122794 01 
Cx —. 640348 —01 . 207820 —-02 . 128917+01 -— 497509 —01 — 142194 . 194136 —01 
Cre —. 531867 . 131177 —01 . 506058 —01 | —. 226239—01 —. 709217 —01 . 870731 —02 
Cu —, 254629 —, 304136 —. 328941 . 875148 —O1 . 124468 —. 101012 —01 
. 194401—01 —. 242335 —01 — 613555 —03 . 532861 —01 — 23418624 —01 
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1 The following convention Is used to Indicate multiplication factors: --O1= 10; -02= x10; —01= 10-1; —02= 10-4; and so forth. 


section 2) represents three times the distance from 180° 
to 120° or from 60° to 0° (corresponding to sections I and 3). 
For the numerical examples, the labor of computation 
necessary to calculate the bending moment and other 
desired quantities may be shortened somewhat through cog- 
nizance that antisymmetrical loading produces an anti- 
symmetrical] moment distribution about the vertica} axis 
and symmetrical loading yields a symmetrical distribution. 
Hence, only half the cylinder at the ring need be considered 
(fig. 4). The procedure used in obtaining the numerical] 
results is illustrated by taking as an example cylinder 2 
subjected to a concentrated radial load at @6=180°. The 
discussion to follow is confined to cylinder 2 so Joaded. 


DIFFERENTIAL EQUATIONS AND SOLUTIONS 
In the top panel, section 1 (fig. 4), the differential cquation is 


aM 
de 


14 d'M 37 d'M 20 dM , 900000 


3 de® 3 d 3 d ' 729 -zzy M=0 





The resulting auxiliary algebraic equation yields the eight 


complex roots (see table JI): 
+ (ea +da) = + (0.915513 +2.1019992) 


Since symmetrical loading is apphed to the structure, the 
moment in this section is given by only the symmetric terms: 


IM ection 1= Cy sinh ¢,6 sin d,6+C, cosh c8 cos d;é-+- 


C, sinh e0 sin da0 + C; cosh & cos da8 (9) 


In the middle panel, section 2, the differential equation is 


d: M d'M a‘ M 60 É á M 


dgs °° ag ae 


—119 = a +900000Af = 0 





The auxiliary algebraic equation has the cight roots (see 
table IT): 


+k, = +5.534415 
+= £7.395633 
+ (a +62) = + (2.181206 + 4.2918822) 
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Elementary theory 


———~—— Cylinder | 


= Cylinder Present theory 





l 
180 120 100 80 60 o 


8, deg 
(a) Ring bending moment. 
(b} Sheet ee flow. 


Fitr 5.—Bending-moment and shear-flow distributions produced by concentrated radial 
load P at §= 190° 
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(a) Ring shear force. 
(b) Ring axial force. 


FIGURE 6.—Ring shear-force and re tears produced by concentrated radial 
at Gm . 
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Elementory theory 


——-—— Cylinder L} Present theory 


Elamentary theory 
Jenss theory 


—— -—— Cylinder 2 





8, deg 8, deg 
: (a) Ring bending moment. r 
ï (bi Sheet shear flow, I Ros sher foice. 


FIGURE 7.—Bending-moment and shear-fow distributions produced by concentrated m°Hent Frayre §.—Ring shear-forca and arial-force distributions produced by concentrated moment 
load Moe at dm 180°, Inad Mfo at 8m 180°. 
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———— Elementary theory 


—————Gylinder >} 
——-——Cylinder 2 Present theory 





§ deg 


(a) Ring bending moment. 
(b) Sheet shear flow. 


Fisvak ¥—Bending-moment ane shear-flow distributions produced by concentrated 
gential load Tat @= 180°, 


The moment in section 2 is then 
hy = Cy’ sinh k8 t+ Cro” cosh k+ Cu’ sinh k+ 
Cy’ cosh k84 Ca sinh af cos b8 + 
C'u sink a8 sin b0+-Ci cosh ag sin 66+ 
Cy, cosh af eos b8 
but since k,@ and k, when @=60° and @=120°, are large 
enough so that their respective hyperbolic sines and cosines 
are of almost equal magnitude, it is better for computational 
purposes to rewrite the moment in section 2 as 
M section a ™ Coet! +- Cpe Fo + Cet + Cpe #2? + 
Cia sinh af cos 6¢+C, sinh ap sin 68+ 
C's cosh af sin b84 Cis cosh al cos b8 (10) 
For the bottom panel, section 3, the differential equation 
is the same as for the top panel and the auxiliary algebraic 
equation has the same roots. Hence, the moment is given by 
NE tiani = Cir sinh ¢,@ cos d+ Oe sinh c,8 gin d0 + 
C,,’ cosh ¢,8 sin dig -+C cosh c18 cos di0 + 
Ca sinh c8 cos d0 + Ca sinh ¢26 sin da0 + 
Ca cosh tof sin d0 +C cosh c.6 cos da? 
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a) Ring shear force. 
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Fiaure 10,—Ring shear-force and axial-force distributions produced by concentrated 
tangential load T at ê= 190°. 


However, for 0=120° and 0=180° the hyperbolic sines and 
cosines, respectively, of ¢,@ are almost identical, and it is 
advisable for computational purposes to rewrite the moment 
in the bottom panel in the form 


M section == C'17€°1* cos d)6+C13¢-t" cos di + 
Ce? sin di8 tC ge 71? sin di8 + 
Ca sink £28 cos df +Ca sinh c0 sin d,6+ 
Ca cosh c0 sin da0 +C., cosh c0 cos d8 (11) 


BOUNDARY CONDITIONS 

Equations (9) to (11) contain twenty unknown constants; 
hence, twenty boundary conditions are needed. No bound- 
ary conditions are found at @=0, since continuity of ali 
forces and displacements is already satisfied as a consequence 
of taking advantage of the symmetry of the structure and 
loading. The boundary conditions to be used in the calcu- 
lations must be found at @=60°, @=120°, and @=180°. 

Eight of the boundary conditions are supplied at @=60°, 
where, in the absence of any concentrated loads, all eight 
expressions in section 1 must equal the corresponding expres- 
sions in section 2. That is, the moment, shear force, dis- 
placements, and so forth must be continuous. For example, 


820 

the first condition requires continuity of the moments or 
M roction 1 — M section a= 0 

which can be written 


LTEC, c0 sh ar cos dr, 


C's sink = sin tery Cs cosh “2% cos hr 


C. sinh ar sin 


kir P a 
Cee M: +C re + is T +Oue F ; 
Cs sinh = eon a teu sinh 22 = T sin ary 
Cs cosh ain! tC cosh 2 = cos sz =0 


The other seven boundary conditions appearing at 0=60° 
are written in a similar manner. The same continuity of 
the forces and displacements must exist at @=120°. Thus 
sixteen boundary conditions are provided, with the four 
remaining conditions to be found at @=180°. 

For an applied radial load, the sum of the shearing forces 
on either side of the cut must equal the applied shearing 
load P. Because of the symmetry of the structure and load- 
ing, half of this load is carried by section 3 and the other 
half of the load is carried by the panel to the left, which 
need not be considered. Furthermore, the shear flow in the 
sheet is 0 at 6=180°, there can be no tangential! or sidewise 
displacement at 6=180°, and there are no angular displace- 
ments caused by a concentrated radial load at this point. 
Hence, the four boundary conditions needed may be summed 
up as follows: 


V ise = 0.5P (12a) 
Qis = 0 (12b) 
Togo ~O (12c) 
dw, vN _ 
ek Be TO (12d) 


Equations (12a) to (12d) can be reduced by proper 
substitution to 


(SF). =0.5PR, 
ic 

(Se OSPR 

a „„=0.5PR, 

(Se go OSPR, 


For example, equation (12a) can be written 
Cne“ (ci cos dir — disin dir)— Cre ™ "le cos dr +d, sin dir) + 
Cie“! (e: sin dir +d, cos dir) — Cne ™ (e sindir — d cos dir) + 
Ca (C: cosh Car COS dst — da sinh Cot sin dr) + 
Ca(C; cosh car sin dert- da sinh car cos dar) 

a(€ sinh car sin dam+ d; cosh cor cos dar) + 
Ca(ca sinh cr cos dzr—d cosh car sin dar) =0.5PR; 
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The solution of the twenty equations given by the boundary 
conditions yields the values of the unknown coefficients Ca, 
C, ... Cy, Cy shown in table III. The bending moment 
in the loaded ring at any angle 8 is then found by using 
equation (9), (10), or (11) with the appropriate coefficients. 
For example, substituting coefficients Ciy to Cy from table 
IIT into equation (11) for 6=180° gives 


Misw= — 0.119054 X 10- PR (2.656444 X 10°) (— 0.372983) + 
(—0.202194 X 10'PR,) (0.376443 X 10-*)(— 0.372983) +- 
(0.342801 X 10-*PR,) (2.656444 X 10°) (0.927839) +- 
(0.216297 X 10'PR,) (0.376443 X 10) (0.927839) + 
(—0.497509 X 10-'PR,) (8.444906) (0.949098) +- 

(— 0.226239 X 10 PR) (8.444906) (0.314983) + 

(0.375148 X 10~1PR,) (8.901257) (0.314983) + 

(0.582861 X 10~'PR:} (8.901257) (0.049098) 
=0.161390PR, 


CONCLUDING REMARKS 


The results obtained from the numerical examples agree 
with those previously obtained for circular cylinders (sce, 
for example, references 2 and 3) in indicating that concen- 
trated loads applied to flexible rings produce stresses in the 
rings and shell considerably different from those computed 
from an enginecring analysis (wherein the ring is treated as 
a free ring supported by the usual elementary torsion and 
bending shears). Ring flexibility is essentially indicated by 


the parameter A (= on) where A is everywhere less than 


about 200 the engineering analysis is adequate, but if A 
exceeds 200 such an analysis is inadequate. 

The main effects of flexibility are to change considerably 
the distribution of stress and the magnitudes of the maximum 
stresses. The change in maximum stresses is indicated in 
the following table, which gives the approximate ratios of 
the absolute values of the maximum stress coefficients for 
cylinder 2 (A +206 near the load but 150,000 some distance 
away) to the maximum stress coefficients obtained from an 
engineering analysis: 





For much larger values of A—that is, for greater flexibility— 
the ratios greater than unity would increase considerably, 
those less than unity would decrease somewhat, and the 
ratios of unity would remain unchanged. 


LANGLEY AERONAUTICAL LABORATORY, 
NATIONAL Apvisory COMMITTEE FOR AERONAUTICS. 
LANGLEY Fye tp. Va., August 7, 1967. 


APPENDIX 
DEVELOPMENT OF GENERAL DIFFERENTIAL EQUATION 


Consider the two-bay panel with uniform geometrical 
properties (fig. 1), loaded in the plane of the middle ring. 
In accordance with the assumption that the shear stress in 
the sheet may vary circumferentially but remains constant 
in the axial direction, an infinitesimal element of the sheet 
is subject to the stresses shown in figure 1. 

Equilibrium in the z-direction requires that 


ot ST tds dz+97 Varið (A1) 
Since 
ri=q 
and 
ot’ =p 
equation (Al) becomes 
Op , OG _ 
oro” (A2) 


Integrating equation (A2) with respect to x gives 
2+ f,(8) 


l 
PEE 


However, since the end ring supports are free to warp out 
of their planes, p=0 at x= L, and therefore 


_7 Of 

Then 
r_n OF 
p=(L— r) = 


The strain relation in the z-direction gives 


u sp 
dr “TEY 
so that 
Ou 


l Ò 
TEAT (A3) 


Integrating equation (A3) with respect to x gives 


u= (Ze 


but since, from symmetry of the structure and loading, 
u=0 at r=0, 


ZN OA E fae) 


f2(#)=0 
so that 
x*\ og 
TEY y (Le 2ds (A4) 
The shear strain can be expressed as 
E _O% , or 
TOE as Oe (A5) 


from which u can be eliminated and v can be found in terms 
of q. Differentiating equation (A4) with respect to ¢ and 


substituting the result into equation (A5) gives 


oe òg 
a GE Et 7 (Ee 2) a8? 
Integrating with respect to x yields 
gq Le r’ 
=i tep (aa) oat 


but since the end ring supports are rigid in their planes, 
v=0 at r= L, and 


L òq 
Je) =- q+ SE Qe? 
so that 
ata 1 o òq 
U —"G I~ SRF (3Lr r — 2L’) Jg 


This displacement in the sheet at the ring, where z=0, 
becomes l 
L? d 
"==; 0+ SRR oe (48) 
The relation between the strain in the sheet in the cireum- 
ferential direction and the deflections at the ring is 


=F z (Sew) 


where by continuity w is the radial deflection of both ring - 
and sheet. Since the strain e in the sheet is the strain at the 
outermost fibers of the ring, 





Thus the radial deflection can be expressed as 


=i hee dr By 





(47) 


Differentiating equation (A6) and substituting the result into 
equation (A7) yields 


(RM, Ldq, D dq 


=-EI Gide | 3IPEF dé (A8) 





The relation between the ring bending moment and the 

radial deflection is given See for example, reference 11) 
oe 

Mp=— fare A oat w) (A9) 


Performing the indicated differentiation of equation (A8) 
and substituting the result into equation (A9) eliminates w 


and gives 
d?* Af, d? 
=p pr (“age Me) ae (ag tgh)+ 


rs d'a +e] 
SEE? do 


Me=— rpa 


(A10) 
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It is now convenient to find the relation between the shear 
flow in the sheet and the moment at the sheet center line, 
as well as to express the moment. at the ring center line in 
terms of the moment at the sheet center line. Figure 3 
shows an infinitesimal segment of the loaded ring and the 
positive directions of the forces. For the first of the three 
equations of equilibrium, taking moments about the origin 
gives 

—dM+RdH+R*h dé+Rm d6+2R'¢ dé=0 


or 
1 (dM R p Rh—Rm) 


=a (d E d (All) 


Summation of the tangential forces (neglecting terms of 
ligher order) gives 


¥de—Rh de—2R¢ de—di =0 


or 
V=% Rh42Rq (A12) 
Summation of the forces in the radial direction gives 
—Hdé+Rkfde—dV=0 
H=— Ss RS (A13) 


Substituting the value of mn from equation (A12) into equa- 
tion (A11) yields the shear force at @ . 


(A14) 


Differentiating equation (A14) and substituting the result 
into equation (A13) gives the value of the normal force at 
the sheet center line at 8 


1 ËM 


H=; Ge oe pe a (A15) 


Differentiating equation (A15) and substituting the result 
into equation (A11) gives the expression for the shear flow 
in the sheet at 0 


map (a tar [Mata (mtae) ew 


The bending moment about the ring center line can be seen 
to be composed of the moment M and the product of the 
axial force H and the eccentricity e or 


Mp= Mf —eH 


e dM 


=M+5, —eRf— (A17) 


Differentiation of equations (A16) and (A17) and sub- 
stitution of the result into equation (A10) leads to the 
desired differential equation 





dE Af deM 6c = dM 
Fee T 3B) Wee +(1- 6.B-b 9 Pr 
2e4 e ENIEM n ee 
3| B— - (2-5ta) | EU aa( $48) arene 
(A18) 
where 


FO=R | a+- 3B) oh a -38 p |+ 


2e ne) Tat 
der 








e| gata) A det 


a (Aft Re) s |+ 
6e% dom 


ae —3 B) Tati- —6B+ -py AET 


| B— “Rt e) | Ta} 


Equation (A18) is the general form of the differential equa- 
tion for the bending moment at the sheet center line for any 
of the panels composing the structure. The bending 
moment in the ring for any such panel is given by equation 


(A17). 
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